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Heisenberg ladders are investigated using the bond-mean-field theory [M. Azzouz, Phys.Rev.B 
48, 6136 (1993)]. The zero inter-ladder coupling energy gap, the uniform spin susceptibility and 
the nuclear magnetic resonance spin-relaxation rate are calculated as a function of temperature 
and magnetic field. For weakly coupled ladders, the energy spectrum vanishes at incommensurate 
wavevectors giving rise to nodes. As a consequence, the spin susceptibility becomes linear at low 
temperature. Our results for the single ladder successfully compare to experiments on SrCu203 and 
(VO)2P2C>7 materials and new predictions concerning the coupling to the magnetic field are made. 
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Among materials that contributed to the increased interest in low-dimensionaLspin systems, one finds ladder 
compounds like vanadyl-pyrophosphate (VO)2P207El or Cu-0 layers of SrCu203.Eftj These can be modeled by the 
Heisenberg Hamiltonian on a ladder geometry. Here, the general Heisenberg model 

H = J± &i,2j ■ &i,2j+l + J'x Si,2j+1 ' Si,2j+2 + J J]] Sjj • Sj+ij — h Sfj (1) 

i,j i,j i,j i,j 

for an array of coupled ladders in a uniform magnetic field h along the z-direction is considered. The quantities 
J, J± and J' ± stand respectively for intrachain, transverse intra-ladder and inter-ladder positive exchange coupling 
constants. It is now well established that ajw finite transverse coupling opens an energy gap, E g , in elementary spin 
excitations for the single ladder (J^_ = 0)in3 As a consequence of this energy gap, thermodynamical quantities like 
the uniform susceptibilty, x', and the nuclear magnetic resonance (NMR) spin-relaxation rate, , show an activated 
behavior at low temper atures,tl T < E g . Recently an important numerical work has been reported by Sandvik et al.u 
who used Quantum Monte Carlo (QMC) technique to calculate the magnetic susceptibility and the NMR relaxation 
rate as a function of temperature for the SrCu203 material. However, comparisons to experiments did not lead to a 
clear conclusion concerning the magnitude of the longitudinal, J, and transverse, Jj_, exchange couplings. The best 
fits to Azuma et alB results they reported for \' an d T^ 1 were obtained using different values, namely J = Jx = 850 
K for x' an d J = J± = 1200 K for 1/T\. It is then interesting to address the question of the relevance of the 
Heisenberg ladder for this system. Frischmuth et alS3 numerically studied (also using QMC method) the temperature 
dependence of the uniform susceptibility for Heisenberg ladders with up to 6 legs. They came to the conclusion that 
the gap rapidly decreases when the number of legs increases (for even number of legs). ___ 

The Hamiltonian (Q) is studied by the bond-mean- field theory (BMFT) explained in detail in Refe-jiBB]. In this 
approach, ([!]) becomes a spinless fermion Hamiltonian 



id j.2j i,2j i,j 



hJ hJ t,3 



using the two-dimepaional (2D) generalization of Jordan- Wigner transformation where the phases $, <&j^ and <J>^ 
are readily derivedO The BMFT treatment consists, first and foremost, to approximate the sum of these phases 
on each plaquette by rr. Then, quartic terms corresponding to Ising interactions are decoupled by introducing the 
bond parameters Q — (ci,jc\ +1 ,■), P = {ci,2jc\ 2j+i) an£ i P' = ( c i,2j+ic| 2^+2) that are calculated by minimising the 
mean-field free energy. The dispersion relation is found to be 

E±(k) = ±^{4Ji sin 2 k x + 2J ±1 J ±2 cos 2 k y + ,]\ x + j\ 2 + 4Ji(Jj_i - J± 2 ) sink x sink y } 1/2 - h (3) 
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with Ji = J(l + 2Q), Jxi = Jj_(l + 2P) and J ±2 = + 2P'), 

In this work, elementary spin excitations are described in a fermionic like picture. Consequently, the number of 
fermions in excitations that contribute to thermodynamical functions is found to be conserved just as in normal 
Fermi systems. This result indicates that only the S z = component of the triplet excited state is relevant to 
thermodynamics, and is.si.milax to what has been already verified both experimentally and theoretically in the case 
of spin-Peierls materials .EiJ£3ii3 For example, this property is reflected by the fact that the magnetic Zeeman term is 
missing in the magnetic field dependence of the zero inter-ladder coupling energy gap E g (T, h) = | Jj_i — J12I/2. The 
simple physical reason behind this effect is that thermal fluctuations cannot create spin configurations with S z = ±1 
because up and down spins are equally distributed. This also remains true at zero temperature because (S z ) = due 
to quantum fluctuations. The energy gap, E g (T,h), of the single ladder (J' ± — 0) reported in Fig. as a function 
of temperature never vanishes, but instead decreases with increasing temperature. In other words, the effect of the 
gap is washed out by thermal fluctuations when the temperature becomes larger than a characteristic value of the 
order of the zero temperature energy gap. ThiSris in perfect agreement with the absence of a phase transition in the 
two real ladder materials we mentioned above.tS Also, from Fig. [j], it is seen that E g (T,h) remains unaffected by 
the magnetic field at low temperature. As we explained before, this result is reminiscent of the fact that only the 
S z — component of the triplet excited state will contribute to thermodynamical quantities, and is a consequence of 
the absence of the magnetic Zeeman term in E g (T, h). The magnetic field dependence is implicitly included in the 
parameters Q, P and P' and vanishes for T — > as can be noticed in Fig. [jj We should mention here that in the single 
laddepJimit (J^_ = 0), the gap at k x = is the same as that at k x = n. This contrasts with exact diagonalization 
resultO where the gap at k x = is twice the gap at k = ir. Our mean-field gap finds however a compromise by being 
half way between these gaps. 

As for inter-ladder coupling, from Eq. (^), one can distinguish three regimes: i) J'^ — where the only transverse 
modes are k y = 0, tt. The third term under the square root in the dispersion relation, that becomes equivalent to 
that of the ladder system, vanishes, ii) J'^ — J±, if Jj_ ~ J then the AF order parameter has to be included.EJ In the 
controversial limit J± -C J, the ground state is either ordered antiferromagneticallyc2l or is a gapless spin liquid state. 
The later is well described by Eq. ([|) where (J±i — J±2)smk x smk y — since J±i — J±2- Hi) J± *C J±, this is an 
interesting regime. The dispersion relation is found to present nodes at incommensurate wavevectors. The third term 
under the square root, which is absent in cases i and ii, is responsible for a such behavior. For J' ± — 0.01 for example, 
the zeros of the spectrum occur at (k x ,k y ) — (3.661212,±7r/2), (k x ,k y ) = (5.763566, ±7r/2) and their images with 
respect to the symmetry centre (tt, tt) as shown in figure |2j. 

Now, we calculate the momentum transfer- and frequency-dependent response function, x(q, ui), by evaluating the 
space and time Fourier transform of the correlation function (T(S Z , ji(t)S z ■ We obtain: 

( n \ 1 Wj /+(k + q)-/+(k) , , /_(k + q)-/+(k) , \ 

= 8^1 { g + (k + q)-£? + (k)-a, + (+ ~* ))6+ + { ^(k + q)- g+ (k)-c + (+ ^ W 

where f± is the Fermi factor and 9± = 1 ± cos(ak — cx-k+q) with otk, satisfying tana*; = 2 Jl ^Tjf^jj^j cofk S1 " ~ > 
is the phase of the coherence factors required in the diagonalization of the mean-field Hamiltonian. The static 
and uniform spin susceptibility, x'(T, h), is obtained by evaluating lim q _ + o.^-+ox(qi ui). This yields the following 
expression: x'(T, h) — & ^ k p=± (l + e l3Epk )~ 2 e~ l3Ep ( k \ The relaxation rate is calculated using Moriya'sEl expression 
= 2j 2 N (gn B )- 2 Tuj- 1 f ^a_|^(q)| 2 x"(q,w) (w -> 0) where x"(q,a;) is the imaginary part of x(q,w) and A(q) is 
the nuclear hyperfine form factor that will be determined by fiting experimental data, /is, g and "fN are respectively, 
the Bohr magneton, the Lande and gyromagnetic factors. In Fig. [I], we display x'(T, h) as a function of temperature 
for several values of the magnetic field for the single ladder limit. From this figure, one can notice that the magnetic 
field manifests itself only at low temperature contrary to the energy gap that is rather affected at high temperature. 
The dominant contribution to x'C^O when T <C E g that is given by cosh(/3/i)e -,3 - E9 (E g is field independent for 
T ~ 0) explains such a behavior. In order to confirm these predictions in experiment, we suggest to carry on a 
measurement of the spin susceptibility as a function of temperature and magnetic field on a ladder material. Note 
finally that a sizable effect of the magnetic field begins at h ~ 0.1E g (T — 0). For weakly coupled ladders, as 
shown in Fig. || for J' ± = 0.01 and J± = J = 1, the energy spectrum is found to show nodes (zero energies) at 
incommensurate wavevectors, (k x ,k y ) = (3.661212, ±7r/2), {k x ,k y ) = (5.763566, ±7r/2) and their images. As a result 
of these nodes, x'{T) ~ T at low temperature. Note that in a work reported by Gopalan et al.j^B the gap vanishes only 
for J' ± / Jj_ > 0.25 at commensurate wavevectors. However, as their mean-field theory favors the opening of an energy 
gap, the question of the right threshold value (wether it is or 0.25) will remain open. An argument in favor of our 
result is given by the fact that nodes occur at incommensurate wavevectors as a consequence of coupling anisotropy. 

Let us compare with experiment. First we would like to address the temperature dependence of the energy gf^i 
In Fig. 0, one notes that the gap decreases with increasing temperature. So, contrary to previous investigations^^ 
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where the zero temperature value of the gap is used for all temperatures, one should take into account the decrease of 
the gap. The data of Johnston et aln concerning the uniform magnetic susceptibility for the (VO)2P2Qf compound 
are fitted in Fig. || using the single ladder limit. These data have already been fitted by Barnes and RieraHl using exact 
diagonalization with J = 7.82 meV and J± = 7.76 meV. We have noted that our calculations for the single ladder 
are more accurate in this case than for SrCu2C>3 (we do not report the fit for SrCu2C>3 because of lack of space). 
Concerning the NMR relaxation rate, we plot in Fig. || both Azuma et al.'s experimental results of the SrCu203 
compound and the best fit obtained using the same coupling constants as for the susceptibility. This suggests, as an 
answer to the problem rose in Dagotto et aVs work,Q that both \' an d V^i can be modeled using the same set of 
parameters of the Heisenberg single ladder. For the magnetic field, a sizable effect in this material is expected for 
fields h ~ Q.\E g (T = 0) <~ 100 T that are not accessible in experiment. One can instead use the (VO)2P20y material 
since the energy gap E g ~ 50 K is much smaller OBil This yields h ~ 10 T. 

The departure from experimental data at high temperature may be attributed to the fact that our spectrum is the 
same for either k± =j-0 or k± = it, and gaps at k — and k — n are equal. As stressed before, the gap at k = is 
twice that at k = 7r.t3 Another point is that since energy gaps decrease with increasing temperature, contributions 
from the continu-am of the spectrum above the gap at k = could become not negligible already for temperatures 
T ~ E g {T — 0).E3 This is obvious for (VO)2P207 where deviations start at around T ~ 50 K. However, it is less 
clear for SrCu203 where deviations are found to occur at T ~ 250 K which is almost half the gap [E g w 468 K). 
From these above remarques, we conclude that our theory can be applied to explain experiments at low temperature 
(T<E g (T = 0)). 

In conclusion, we have used the bond-mean-field theory to investigate AF Heisenberg ladder spin systems. The 
uniform spin susceptibility and the NMR spin-relaxation rate are calculated as a function of temperature and uniform 
magnetic field. For a single ladder, the energy gap remains unaffected by the magnetic field at low temperature 
in agreement with the fact that the relevant elementary spin excitations to thermodynamics (as seen here for x' 
and Tf 1 ) conserve the total S z = component. Our results successfully compare to experiments on SrCu203 and 
(VO)2P207 compounds for T < E g . Several predictions concerning the magnetic field effect that can be verified in 
simple experiments are also made. Inter-ladder coupling is found to drastically change the behavior of the uniform 
susceptibility, at very low temperature. For J' ± = O.OlJi for example, \' becomes linear in T as a consequence 
of the appearance of nodes in the spectrum at incommensurate wavevectors, (k x , k y ) — (3.661212, ±7r/2), (k Xl k y ) 
(5.763566, ±7r/2) and their images. It is worthy of note to mention, in the end of this work, that Gopalan et al$3 
suggested that the inter-ladder coupling in SrCu203 is ferromagnetic. They based their argument on the fact that two 
Cu ions on adjacent ladders are connected through an O site by 90° bonds. A calculation in this case of inter-ladder 
coupling will be reported in the near future. 
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FIG. 1. (a) The susceptibility is drawn for J' ± — as a function of temperature for several values of h. In the inset, E g is 
drawn as a function of temperature for h = 0. 



FIG. 2. The dispersion relation for weakly coupled ladders is ploted as a function of k x and k y ; J' ± — 0.01J±. 



FIG. 3. In (a), the fit to Johnston et al. susceptibility data is reported for (VO)2P207. The fit to Dagotto et al.'s experimental 
results for \ jT\ is shown in (b) in the case of SrCuOs. Data are drawn in crosses. 
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